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I ' Abstract 

(— I \ The significance of NNLO (3-loop) QCD contributions to tlie flavor non-singlet 

sector of F!^^ and F^ has been studied as compared to uncertainties (different 
factorization schemes, higher twist and QED contributions) of standard NLO (and 
r> \ LO) QCD analyses. The latter effects turn out to be comparable in size to the 

C^ ' NNLO contributions. Therefore the minute NNLO effects are unobservable with 

presently available (precision) data on non-singlet structure functions. 



1 Introduction 

In a recent publication ,lj a next-to-next-to-leading order (NNLO) QCD analysis of 
F2^{x, Q"^) and F2'^{x, Q"^) in the flavor non-singlet sector was presented. Here our purpose 
is to study the significance of the NNLO contribution as compared to other, possibly 
important, contributions such as redundant terms in the NLO analysis (which arise when 
the NLO evolved parton distributions are multiplied by the coefficient function), higher 
twist and QED contributions to NLO, and effects due to choosing different factorization 
schemes. In Sec. 2 we present the relevant theoretical expressions required for our analysis, 
and Sec. 3 contains the quantitative results. Our conclusions are summarized in Sec. 4. 

2 Theoretical Formalism 

The non-singlet (NS) parts of the structure functions ^2^' (x, Q^) for x > 0.3, where 
valence quark dominance is adopted, are, at LO, given by 

F2 =-XUv + -Xay = — ^ ^NS,8 + 7^ ^ 9'ns,3 \^) 
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F2'' = -^ xiu, + 4) = Y^x g+s,8 (2) 

where d = {p + n)/2 and g^ss = '^v — dy. For x < 0.3 one analyzes the genuine NS 
combination 

F|-" ^ 2(Ff - Ff ) = i xiuy - ci.) + ^ xiu -d)^lx q^,^, (3) 

where now g^g ^ = u^ — dy + 2(u — d) since sea quarks cannot be neglected for x smaller 
than about 0.3. For definiteness we adopt foid — u the choice |21Q] 

x{d-u){x, Ql) = 1.195x^-2^(1 - a;)^-^(l + 14.05x - 45.52a;2) (4) 

at Qo = 4 GeV^ which gives a good description of the Drell-Yan dimuon production data 
[S], but plays a marginal role in our analysis. 



At NLO(MS) the n-th Mellin moments of the above NS combinations of valence parton 
distributions, for brevity denoted by v^, symbohcally evolve according to the well known 
expression (see, e.g. [HE]) 

v+{Q') = {1 - (a - ao)Ri} (a/ao)-<s//3o ^+(g2) ^5) 

where R^ = PSi^/Po - (/5i//3o)^ii°s and a = a{Q^) = a,(Q2)/4vr with ao = a{Ql). The 
moments of the above NS structure functions F^^ are then given by 



.(1) 



F,^''{Q')=[l + aC'^^v^Q') (6) 

and this expression is commonly compared with experiment. Inserting f^(Q^) from (5) 
into this equation one observes a redundant O^a^) contribution, i.e. —a{a — ao)C2 j(fgi?i, 
which in fact belongs to a NNLO analysis and is assumed to be small at NLO. If, however, 
one chooses to work to a NNLO (3-loop) accuracy, such a redundancy at NLO might 
become significant as compared to the full NNLO contribution. This will be investigated 
quantitatively below. Similarly, the choice of a factorization scheme, other than the MS 
scheme used thus far, might imply larger differences than additional NNLO contributions 
in the MS scheme. For example, in the deep inelastic scattering (DIS) factorization scheme 
[HI HI the Wilson coefficient in (6) is absorbed into the parton distributions, i.e. into their 
evolutions in (5) with Ri -^ i?f ^^ = Ri - C^%s '■ 

and, instead of (6), F^^ = i;^jg at NLO. 

Taking into acccount QCD 3-loop NNLO af effects, one also has to consider QED LO 
a-contributions which are of comparable size |7|. The latter can easily be implemented 
by changing the n-th moments of the NLO valence input distribution QviQl), Qv = u^, d^, 
in (5) according to 
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QviQo) ^ {a/aoY^ "^ exp 
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with a ~ ^ and P^ = (e2/CF)Piis where Cp = 4/3. 

At NNLO(MS) the evolution of v+{Q'^) in (5) generahzes to 
v+{Q^) = {l-{a-ao)Ri-^{a^ -al){R2-Rl)-ao{a-ao)Rl}{a/ao)-<^^^' v+iQl) (9) 
with i?2 = PJ,T/Po - {Pi/Po)Ri - {P2/P^)P}S and (6) becomes 



^(1) _L ^2ry{2)+ 



F.^'iQ') = [1 + aC^:^s + ^'C'2%\ ^^m ■ (10) 

('21+ 

Convenient expressions for the relevant 2-loop Wilson coefficient C2 ^s ^^^ be found 
in |5j (eq. (A. 2)) and for the 3-loop sphtting function P^^ in ^ (eq. (4.22), which 
can be easily Mellin-transformed using 0). The strong coupling now evolves according 
to da/dlnQ"^ = -ELo/^^«^"^^ ^^^re /5o = 11 - 2//3, pi = 102 - 38//3 and p2 = 
2857/2 — 5033//18 + 325/^/54, which refers to the MS renormalization scheme, and / 
denotes the number of active flavors. Here the redundant (9(a^,a^) contributions in (10) 
turn out to be marginal and do not influence our (flt) results. Furthermore, the running 
coupling a(Q^) is appropriately matched at Q = m^ = 4.5 GeV and Q = rrit = 175 GeV. 
To obtain F|^^ in the DIS factorization scheme at NNLO, the Wilson coefficient functions 
in (10) have to be absorbed into the parton distributions, i.e. into their evolutions in (9) 
with R2 ^ R2 = R2 — 2C2 Ns + (^2 Ns) • 

1 
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,m = {l-{a-ao)R?''-;-{a'-al)[Rr-{Rrr] 






where i?f^^ is as in the NLO-DIS expression (7) and, instead of (10), we now have 
^2^3 _ y+ ^ ^^ NNLO. 

Since flavor NS structure functions are mainly related to the medium and large x- 
region, the relevant kinematic nucleon target mass (TM) corrections are always taken 



into account according to [TU] 

^2,TM\'^^Q ) = I -^ ^2,TM\-^^Q ) dX = 

Jo 

^WV jKn-2)\ {n + 2j){n + 2j-l)^ \\Q'J j^ ' 

where higher powers than {m^/Q'^Y ^^^ neghgible for the relevant x < 0.8 region, as can 
straightforwardly be shown by comparing (12) with the well known exact expression in 
Bjorken-x space jlOj. 

Despite the kinematic cuts (Q^ > 4 GeV^, VT^ = (i - 1) g^ + ^2^ > ig GeV^) used 
for our analysis, we also take into account higher twist (HT) corrections to Ff'^ via 



v^{x,Q^ 



m% 



l + ^hix) 



v''{x,Q') (13) 



Q 

in order to learn whether nonperturbative effects may still contaminate our perturbative 
analysis. Here we adopt the ansatz ^I] 

h{x) = a(— cj . (14) 

Notice that the input valence parton distributions v~^{x,Ql) at LO [v^{Q'^) = 
(a/ao)~-^Ns//3o v+{Ql)], NLO and NNLO in the MS as well as in the DIS scheme in eqs. (5)- 
(11) and (13) can and will be different in general. 

Finally, Fermi motion and nuclear effects in the deuteron are strongly model dependent 
and will therefore not be considered here. They were, however, taken into account in 
mElEl using the specific models cited there. Comparing these results with our valence 
distributions obtained and to be discussed below, as well as with other results where such 
effects have not been taken into account (e.g. [2]), shows that these effects do not change 
the quality of the QCD fits. 



3 Quantitative results 

In the present analysis we used the proton and deuteron data of BCDMS fT^, NMC 
[T^ and SLAC 15 , as well as the proton data of HI TBj and ZEUS |12| in the relevant 
x-regions discussed above which amount to 480 data points. The valence distributions 
have been parametrized at the input scale Ql = 4 GeV^ as 

xu^{x,Ql) = NuX^^il-xY^il + AuX^^+Bux) (15) 

xd,{x,Ql) = Ndx'''{l-xp{l + Adx''' + Bdx) (16) 

with the normalizations A^^ and A^^ being fixed by J^ u^dx = 2 and J^ d^dx = 1, respec- 
tively. The LO, NLO and NNLO fit results without HT contributions are summarized 
in Table 1. The standard NLO and NNLO fits refer to the MS scheme according to 
eqs. (5), (6) and (9), (10), respectively. Our fit results for NNLO are compared in Fig. 1 
with the data used. Except perhaps in LO, we obtained equally good and acceptable fits 
(x^/dof) in each perturbative order and scenario. As has been already noted previously 
Pl El QHl Onj, a NNLO analysis in general results in a slightly smaller as(m|) than in 
NLO. This is due to the fact that the higher the perturbative order the faster as(Q^) in- 
creases as Q^ decreases. In order to compensate for this increase, a NNLO fit is expected 
to result in a smaller value for asi'm'^z) than a NLO fit. Notice that the values of as{Tn?z) 
obtained in the usual perturbative NLO and NNLO fits in Table 1 are comparable to the 
ones in ^1^1^]. Repeating the NLO and NNLO fits in the DIS factorization scheme 
improves only marginally the global MS fits (x^)) and the QED 0{a) contributions leave 
the original NLO (MS) results practically unchanged as evident from Table 1. 

On the other hand, the inclusion of higher twist contributions sizeably improves the 
fits, i.e., the value of x^/dof as can be seen from Table 2 {as{m?z) is reduced as expected 
since the HT term takes care already of some of the Q^-dependence of the data). In order 
to illustrate the relative significance of the (model dependent) HT corrections in (13), we 
have performed a fit for Q^ > A GeV^ and one for Q^ > 10 GeV^, denoted by HT(IO) in 



Table 2. Clearly, HT effects become less important when the lower cut of Q^ is increased 
and the value of x^ increases, eventually approaching the larger values obtained by purely 
perturbative fits. Nevertheless it is remarkable that the fits for Q^ > 4 GeV^ and Q^ > 10 
GeV^ are not significantly different as will be illustrated below (Fig. 2). For illustration 
we also show in Fig. 1 the NNLO results with HT effects included for Q^ > 4 GeV^. The 
results for the Q^ > 10 GeV^ cut are not shown, since they are very similar to the ones 
for the Q^ > 4: GeV^ cut (dashed curves). 

The actual relative size of our results can best be seen by comparing the various fit 
results with the pure QCD NLO fit, i.e. by considering the following ratios, depicted in 
Fig. 2, which are defined as follows. The effect of NNLO contributions can be visualized 
via rNNLO = -^2^nnlo/-^2^nlo with the nominal NLO structure function given by (6) and 
the NNLO one by (10). Similarly, Tlo requires the usual -F2^lo ^^ given by (5) and (6) 
with Ri = and C2 = 0. Furthermore r^^Q = -F2^nlo /-^2^nlo illustrates the effects of 
choosing the DIS factorization scheme instead of the MS scheme, with the DIS structure 
function given by (7). Similarly, the definition of t^JPlq employs the DIS structure func- 
tion in NNLO given by (11). Including the QED 0{a) contributions to -^2^^lo according 
to (8) results in tJ^lq = -F^nlo /-^2Nlo- Figure 2(a) demonstrates that ambiguities of 
standard NLO MS analyses or additional QED contributions are comparable in size to 
NNLO contributions in the relevant medium and large x-region. Only in the smaller x- 
region around x = 0.2, the NNLO results are about 1% larger than the NLO ones. This 
difference, however, disappears in the DIS factorization scheme where r^J^^Q and t^^q are 
comparable and small (about 0.5%). We therefore conclude that the DIS scheme guar- 
antees a better perturbative convergence than the commonly used MS scheme, except in 
the very large x-region (where nonperturbative contributions are uncontrollable anyway). 
Unfortunately such minute effects are not testable with presently available precision data 
for non-singlet structure functions which have a typical uncertainty of about 10% in the 
small and large x-region. 



The redundant O(a^) contribution to a NLO analysis, as discussed after (6), turns 
out to be marginal: repeating the NLO fit with the redundant term removed from -F2^nlo 
in (6), 

one obtains |rj|^ ~ M = I-^2^nlo/-^2^nlo ~ 1| ~ 0.001 for all relevant values of Q^. 

Finally, the relevance of HT effects is illustrated in Fig. 2(b) where the ratios rf^^ = 
^2i /-^2^NLO f*-*^ ^^^ different perturbative orders are shown, with the HT corrections 
incorporated according to (13). In general, for each perturbative order the fit results are 
stable with respect to different choices for the lower bound on Q^ in (13), i.e., Q^ > 4 GeV^ 
and 10 GeV^ denoted by HT and HT(IO), respectively. In particular the more relevant 
NLO and NNLO fit results are very similar. Despite the fact that these results are about 
twice as large as the ones without HT contributions in Fig. 2(a) throughout the whole 
x-region considered, they are still much smaller than present experimental uncertainties. 

4 Conclusions 

The significance of NNLO QCD contributions to the fiavor non-singlet sector of 
F2^' (x, Q"^) has been studied as compared to uncertainties of standard NLO (and LO) 
analyses. NNLO corrections slightly improve the fits to presently available data and imply 
a better perturbative convergence in the DIS factorization scheme than in the commonly 
used MS scheme. However, ambiguities of NLO fits such as the choice of a particular 
factorization scheme (MS vs. DIS) and possible higher twist effects as well as QED 0{a) 
contributions turn out to be comparable in size to NNLO (3-loop) contributions. In par- 
ticular, nonperturbative higher twist effects play an important role in obtaining optimal 
fits (minimal x^) which turn out to be rather stable with respect to different choices of 
the lower bound on Q^. Their contribution is about twice as large as purely perturbative 
uncertainties which are typically less than about 1%. We therefore conclude that the 



rather minute NNLO QCD effects in the ffavor non-singlet sector are not observable with 
present precision data for flavor non-singlet structure functions which have sizeably larger 
errors. 



This work has been supported in part by the 'Bundesministerium fiir Bildung und 
Forschung', Berlin/Bonn. 
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Table 1: Parameter values of the QCD fits in various perturbative orders based on all 
non-singlet data for Fg^^' (x, Q"^). The parameters of the input valence distributions refer 
to (15) and (16). The analysis for the DIS factorization scheme is based on (7) for NLO 
and on (11) for NNLO. The QED contribution to NLO is taken into account according 
to (8). 





Parameter 




xVdof 


as{ml) 




(lu 


bu 


Cu 




(A(4)/MeV) 




An 


Bu 


Nu 








a-d 


bd 


Cd 








Ad 


Bd 


Nd 






LO 


0.574 


3.290 


0.823 


0.98 


0.128 




9.258 


-7.164 


1.790 




(196.0) 




0.600 


4.952 


0.100 








-0.330 


3.726 


1.742 






NLO 


0.600 


3.364 


0.667 


0.93 


0.112 




9.900 


-8.504 


1.521 




(222.6) 




0.600 


5.163 


0.181 








3.294 


8.303 


0.542 






NLO DIS 


0.600 


3.004 


0.855 


0.91 


0.113 




9.801 


-9.010 


2.101 




(228.3) 




0.581 


4.797 


0.878 








-0.030 


4.458 


1.264 






NLO QED 


0.600 


3.361 


0.667 


0.93 


0.112 




9.883 


-8.496 


1.523 




(217.1) 




0.599 


5.161 


0.186 








3.223 


8.037 


0.555 






NNLO 


0.600 


3.571 


0.599 


0.89 


0.111 




9.654 


-6.810 


1.309 




(177.2) 




0.600 


5.209 


0.323 








4.060 


4.870 


0.657 






NNLO DIS 


0.587 


2.727 


0.825 


0.89 


0.112 




9.644 


-9.471 


1.941 




(187.2) 




0.600 


4.787 


0.868 








-2.004 


6.046 


1.421 







Table 2: As in Table 1 but including HT contributions as well according to (13) with the 
parameters (a, b, c) referring to (14). Fits using a lower bound Q^ > 10 GeV^ in (13) are 
denoted by HT(IO), whereas HT refers to Q'^ > 4 GeV^ as stated before (13). NLO and 
NNLO always refer to the MS factorization scheme. 







Parameter 


xVdof 


as{m%) 




a« 


bu 


Cu 




(AW/MeV) 




Au 


Bu 


Nu 








ad 


bd 


Cd 








Ad 


Bd 


Nd 








a 


b 


c 






LO HT 


0.600 


3.218 


0.397 


0.82 


0.118 




6.926 


-4.208 


1.021 




(121.8) 




0.582 


4.882 


0.712 








7.205 


-2.878 


0.956 








2.180 


0.941 


0.577 






LO HT (10) 


0.600 


3.167 


0.422 


0.88 


0.126 




8.286 


-5.214 


0.942 




(176.6) 




0.597 


4.985 


0.750 








6.163 


-1.061 


1.052 








2.105 


1.000 


0.749 






NLOHT 


0.600 


3.368 


0.508 


0.83 


0.104 




9.897 


-8.215 


1.104 




(136.6) 




0.600 


5.383 


0.222 








4.905 


9.563 


0.467 








2.198 


1.468 


0.303 






NLO HT (10) 


0.600 


3.429 


0.522 


0.90 


0.106 




9.900 


-7.786 


1.126 




(155.8) 




0.598 


5.507 


0.190 








2.383 


9.824 


0.642 








2.313 


1.740 


0.284 






NNLO HT 


0.600 


3.571 


0.519 


0.83 


0.103 




9.900 


-7.086 


1.116 




(109.5) 




0.566 


5.425 


0.302 








5.816 


9.475 


0.441 








1.695 


1.104 


0.432 






NNLO HT (10) 


0.600 


3.622 


0.518 


0.91 


0.105 




9.898 


-6.491 


1.099 




(124.9) 




0.599 


5.594 


0.235 








2.260 


9.119 


0.716 








1.989 


1.801 


0.252 







3.0 



2.5 



2.0 



1.5 



1.0 



0.5 



0.0 



J^l x=0.0125(+2.6) 
+^ 0.0175 (+2.4) 
^-H-M 0-025 (+2.2) 
■t ^ j , ) , j I 0.035 (+2.0) 
n ttj, ), M t 0.05 (+1.8) 

0.07 (+1.6) 
t M t I 0.09 (+1.4) 

m t*^ I t 0.1 (+1.2) 

+ t t t n M ' ' 0.11 (+1.0) 

M I, t tt|tj('"|.t4* t 0.18 (+0.6) 
m t | I, ^ ♦^ ♦y*<'> - t'*** 0.225 (+0.4) 

0.275 (+0.2) 
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Figure 1: Comparison of our NNLO fits with all presently available flavor non-singlet data 
[T^ im n^l [T^ [T7j used for our analysis. The higher twist (HT) contribution is taken 
into account according to (13) and (14). The NLO fits are very similar and practically 
indistinguishable from the ones shown. So is the NNLO HT(IO) fit resulting from the cut 
Q^ > 10 GeV^. The inset shows our NNLO input valence distributions at Ql = 4GeV^. 
The scales on the left ordinate refer only to Ff" where for each fixed value of x we have 
added the constant in brackets to -F|'~"'. The scales on the right ordinate refer to F2 and 
to F2. The data sets are shown with their normalization factors in parentheses (first entry 
refers to Ff, second entry to F^) as obtained in the fit. The ZEUS data for F|' have been 
shifted to the right by 5 % in order to make their error bars distinguishable from the ones 
of the HI data. 
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Figure 2: The size of perturbative (a) and nonperturbative (b) uncertainties encountered 
in various perturbative a^-orders of QCD analyses relative to our nominal NLO analysis 
of -^2 NLO (-^5 Q"^) which always appears in the denominator of the ratios r as defined and 
discussed in the text. The HT contributions to the fits shown in (b) refer to the cut Q^ > 4 
GeV^, whereas HT(IO) refers to a cut Q^ > 10 GeV^. The typical relative experimental 
accuracy is illustrated at a; = 0.4 by the vertical bar (±1 %). At larger and smaller values 
of X the experimental error increases (e.g., the uncertainty is about ±2.5% at a; = 0.55, 
and ±10 % at X = 0.18). All results are shown for Q^ = 40 GeV^, but the agreement with 
data at g2 = 4 GeV^ and Q^ = 100 GeV^ is similar. 



